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, Abstract 

[ — , We study periodic arrays of non-Abelian vortices in an SU (N) x U(l) gauge theory with Nf flavors of 

fundamental matter multiplets. We carefully discuss the corresponding twisted boundary conditions on 
the torus and propose an ansatz to solve the first order Bogomolnyi equations which we find by looking 

i-Sh , to a bound of the energy. We solve the equations numerically and construct explicit vortex solutions. 

^ ■ 1 Introduction 

Vortices -string-like configurations with quantized magnetic flux- have been introduced in the context of 
J> ■ High Energy Physics more than 30 years ago in the pioneering work of Nielsen and Olesen [1] on the Abelian 

Higgs model. Very soon it was shown that in this model the equations of motion can be reduced to first 

00 i 

, order Bogomolnyi-Prasad-Sommerfield (BPS) or self-dual equations [2], [3], explicit numerical solutions were 

04 | obtained and the connection with supersymmetry was signaled [3J . 

The formal proof of existence of a general n- vortex configurations for the BPS equations in R 2 (n being the 
vorticity or number of flux-quanta) was later given in [4] . Although the existence of vortex solutions for the 
BPS equations on other (Kahler) manifolds was shown in [5], it was only very recently that Gonzalez Arroyo 
and Ramos [3] presented a simple method to construct explicit numerical solutions on the bi-dimensional 
torus T 2 , the simplest compact manifold. The case of the torus is particularly interesting as it corresponds 
to the study of periodic field configurations, leading in this case to vortex-lattice arrays (a situation that 
most often arise in condensed matter problems). The torus also provides the most natural (long-distance) 
regularization of R 2 , so even for the cases in which the problem of interest is set in the plane, it is important 
to have a well controlled method that allows us to study the asymptotic infinite area case. For instance, the 
torus provides a natural set up to study numerically non-cylindrically-symmetric multi vortex configurations 
(the case of n superimposed vortices or "giant" vortex can be more easily treated on a disc) . 

Self-dual equations are of course simpler to study than the Euler-Lagrange equations but, more impor- 
tantly, it is their connection to supersymmetry and their relevance to the understanding of non-perturbative 
phenomena in field theories and string theory that have triggered so many investigation in the last years (for 
recent reviews on the subject see [7], [8]) . 

It was shown long ago that non-Abelian gauge models also have vortex solutions [5] (see also references 
in [7]-[S]). More recently much attention has been devoted to the analysis of certain M = 2 supersymmetric 
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theories where non-Abelian vortices have been presented and used to construct low-energy effective actions 
related to string dynamics [10]- [16]. With this in mind, local and semi-local vortex solutions have been 
studied in models in which Yang-Mills fields are coupled to matter fields with different numbers of colors 
and flavors (see [T7] for a review). 

In this work we construct non-Abelian vortices in a gauge theory defined on the torus. Motivated by the 
great interest that the type of theories considered in [TU]-[I1] has received recently, we consider a model with 
gauge group SU(N) x U(l) with Nf flavors of fundamental matter multiplets. After defining the model in 
section 2, we discuss in section 3 the appropriate boundary conditions on the torus which will lead to the 
non-trivial magnetic flux associated to the vortices. In section 4 we obtain a bound for the energy per unit 
length of the static vortices (the string tension) and from it we derive the first order Bogomolnyi equations 
(which correspond to a minimum of the energy). We present vortex solutions to these equations in section 5 
for the case N = Nf = 2 and then extend the analysis to general N in section 6. We summarize and discuss 
our results in section 7. 



2 The model 

We consider the bosonic sector of an M = 2 supersymmetric Yang-Mills-Higgs with SU(N) x U(l) gauge 
group and Nf flavors of fundamental matter multiplets, described by the action 13J 

S = Jd 4 x ±^F™F m ^ + ^F^F ^ + (D^^D^f) 



4<7' 

%r I 4>U m 4> f I 2 ' 



(i) 



In the present case, coordinates x\,X2 are defined on the two-torus T 2 of size Li x L2. The scalar matter 
fields are denoted by <j) with / the flavor index, / = 1, . . . ,Nf, n, v = 0, 1,2,3 are Lorentz indices and 
m = 1, N 2 — 1 is an internal space index of SU (N). We indicate with a superindex the U(l) components 
of the gauge fields. The SU (N) generators t rn are taken to be anti-hermitian with the following normalization 

Tr {t m t n ) = --S mn , {t m )) (< m )f = —SiS* + — Sffi , [t n \t n ] = f mnl t l . (2) 

Here i, j — 1,2, . . ., N are indices in the fundamental representation of SU(N). The covariant derivative 
reads 

D^f = 0^ - A™t m <t>f - A^V (3) 

and we take the U(l) generator t° = (i/\ // 2N)I. Field strengths associated with gauge fields — A™t m 
and A° are defined as 

F MV = d^A„-d„A^-[A^A„} = ^L=F°J + F™t m 



'2N 

F° - r) A — r) A 

F™ = d»A™ - d v A™ - r nl AIA\ . (4) 

The two last terms in (fTJ) are responsible for gauge symmetry braking. The last one, containing the 
Fayet-Iliopoulos parameter £, forces to develop a vacuum expectation value, while the last but one forces 
the VEV to be diagonal. We shall be interested in the case Nf > N since (when Nf < N there is spontaneous 
supersymmetry breaking [7]). 
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Up to gauge transformations, the minimum of the potential, which we call <E»o 5 can be written as an 
N x Nf matrix, 

(l o ... o\ 



(5) 



\0 1 ... 0/ 

Here rows correspond to the N colors and columns to the Nf flavors. Such a minimum breaks the symmetry 
U(l) x SU(N) x SU(Nj) down to a global remanent symmetry SU(N)c+f, 

£7(1) x SU(N) x SU(Nf) -> SU(N) C +F ■ (6) 

The global SU(N)c+f transformations leaving the vacuum invariant act on 4>o according to 

®o = U- 1 <S> ^j , (7) 
where U is a global SU(N) matrix, and is the null (Nf — N) x N matrix. 

3 Boundary conditions 

Since we are working in the two-torus (xi, x 2 ) C T 2 , the gauge fields and the matter multiplets <fif must 
obey periodic boundary conditions on T 2 modulo gauge transformations 

A^xi + U) = UiA^Uf 1 + dpUiUf 1 (8) 
tfSfa+Li) = Uitf , (9) 

where f7, G x SU(N) (i — 1, 2) are the transition functions. Consistency of equation ((SJ for the gauge 
fields leads to the following relation for the transition functions 

U 2 (xi +Li,x 2 )[/i(xi,x 2 ) = Ui(xi,x 2 + L 2 )U 2 (xi,x 2 ) • (10) 

Here ft is an element of the Zm x t/(l) center of U(N) [TS] which labels inequivalent topological sectors. 
Its presence is due to the fact that transition functions Ui can be defined in (jSJ) modulo an element of the 
U(N) center. However, in the present model there are matter fields in the fundamental representation and 
one has also to check consistency of equation (J9]). In this case one finds 

U 2 (xi +Li,x 2 )[/i(a;i,X2) = Ui(x 1 ,x 2 + L 2 )U 2 (xi,x 2 ) , (11) 

so that SI in (jTCTjl should be taken as the unity matrix, fl = I. 

In order to find the solution to equation (111)) it will be convenient to construct the following linear 
combination of elements in the Cartan subalgebra of U(N) 

T ° = diag ^ 1 ' " ' ' ^ 

n = . * diag(-(JV-l),l,...,l) 

y/2N(N-l) 



t n = = diag(l,...,l,-(JV- 1)) ■ (12) 

^2N{N- 1) BV V ^ V ; 



3 



We then define a g-elementary transition function solution (x\, x 2 ) as that which is generated by tq 
(which is proportional to the identity) and t 9 , 



U^ J (xx,x 2 ) = Hexp( — — 



U^( Xl ,x 2 ) = nexp(^g) , (13) 



with 



, 2 

70 = -2tt 



„ 2(N-1) 

7l = Vi^0,g. (14) 

Using this transition function, we shall be able to construct elementary vortex solutions in the torus which 
in the L\,L 2 — * 00 limit reduce to the elementary vortex solutions in [llj.|14j. Such elementary vortices 
have a quantum of magnetic flux. 

Vortices with higher units of magnetic flux are constructed by considering transition functions which are 
products of elementary transition functions. These general transition functions are connected with those 
introduced by 't Hooft 18J, 19 . In fact, a general transition function with n units of magnetic flux reads 

7o = -2wJno 



lq = 2 *^ N N ~ n i > q = l,---,N, (15) 

where 

N 

n o = Yl n q . (16) 

q=l 

With the transition functions that we have defined, conditions ® read 

A 1 (xi+Li,x 2 ) = Ai(xi,x 2 ) 

A 1 (x 1 ,x 2 + L 2 ) = A 1 (xi ,x 2 )- — — 7 ro r m 

^ m 

A 2 (x 1 +L 1 ,x 2 ) = A 2 (xi,x 2 ) + — — y^7 m r m 

IL 2 c — ' 

ra 

A 2 (x 1 ,x 2 + L 2 ) = A 2 (xi,x 2 ) , (17) 
and can be written as a sum of a periodic function Ai in T 2 plus a known function 

Ai(xt,x 2 ) = Ai(x 1 ,x 2 ) + 777-7- ^ 7m T m e tj Xj , (18) 

IL\L 2 — 

m 

so that the field strength takes the form 

F - = F ^-Vt t 

F l3 = diAj-djli- [li.ij] . (19) 
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The magnetic flux can be defined from (fTO]) as 

® mag = - iTy J T2 dTF ^ = = 277 n °- ( 2 °) 

In analogy with what we did for the gauge field, we shall propose an ansatz for matter fields factoring 
out in each multiplet a particular N x N matrix \ — (Xij( x ii x 2)) satisfying the twisted boundary condi- 
tions, times a scalar multiplet carrying both color and flavor indices, A((xi,X2) obeying periodic boundary 
conditions 

<p f (x 1 ,x 2 ) = x(xi,x 2 )A f (xi,x 2 ) . (21) 
Given conditions (O, x( x ii x 2) should satisfy 



\(xi+L 1 ,x 2 ) = Y[cxp ( --^Y^-^J X(xi,x 2 ) 



X(xi,x 2 +L 2 ) = JJexpj^gJx^,^). (22) 



A solution to these equations is given by 

N 



X{xi,x 2 ) = [I exp I — —\Q m (x 1 ,x 2 ) (23) 

m— ^ 12 / 

with 

9"0n, x 2 ) = n e 3 ( i ^ mTm {X1 + ix T 2 + a - } |_:m™M . (24) 

^ V 2n m L\ 2-Kn m LiJ 

Here, n m is the charge associated to the m-th direction of the Cartan subalgebra as defined in (p~2|) . and 

3 (z | r) =£V 7rTp+2fc , (25) 
is the Riemann Theta function. The complex coefficients a™ satisfy the conditions 

n=l 

and determine the position of the vortices. 

In the next sections, A* (x\,x 2 ) in (|21j) together with Ai(a;i, #2) will be determined from the equations 
of motion. 



4 Bogomolnyi equations 



We are interested in infinitely long (£3— independent) static configurations which extremize the tension 
(energy per unit length), 



T 



dT 



_L jpm pij j \_ pO pij 

, 2 ± ij ^ m ~ ^2 V ' 



*+ 1 Di<j) f i 2 +^ 1 d>)t^ i 2 +^ (<^/ _ 



(27) 



where dT = dx 1 dx 2 is the integration measure over the torus T 2 . 

Since action ([1]) is the purely bosonic part of an J\f = 2 supersymmetric action, coupling constants and 
the form of the potential are automatically adjusted [20] so that Bogomolnyi completion can be performed 
[2]. Indeed, using the relation 



DnjS | 2 H (D 1 ±iD 2 ) ( j>f | 2 T * \ 4>\ Fij <\>U l3 + 



td 



(28) 
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with "td" a total derivative term, we can write the energy per unit length as 

2 / ., \ 2 



T= dT 
Jt 2 



(29) 



Using equation (|2H)l we find that the tension is bounded by 

T > 2tt I n | £ > ™o£2 (30) 
with no defined in (|16p . This bound is saturated whenever the following Bogomolnyi equations hold 

F% = ±ig 2 ^ f t m ^e ij (31) 

= T ^(^ _iVe ) e « (32) 

n,<y = i ", / /^o- ! . (33) 

As we already mentioned, solutions to these equations will also satisfy the second order Euler-Lagrange 
equations of motion. For dcfinitcncss we shall choose the upper sign in these equations (the other choice can 
be handled analogously). 

We shall look for solutions to the equations ([3"T]) - (f3"3")) subject to the boundary conditions discussed in 
the previous section. We start from eq. ((33|) . 



{Di+iD 2 )4> 1 =0 , (34) 

and write <fr(x\,X2) in the form (|21[) 

cj> s {xi,x 2 ) = x{xi,x 2 )A f (xi,x 2 ) ■ (35) 

We have already found the explicit form for \ (see eq.(|23p). which was obtained by fulfilling the required 
boundary conditions. We shall now determine A-^ so that $ in (f35|) satisfies Bogomolnyi equation (|34| . To 
this end it is convenient to write 

Af(x 1 ,x 2 )=M(x 1 ,x 2 )V f . (36) 

Here V? is a constant multiplet carrying both color and flavor indices, while M[x\,x 2 ) € SU(N) x U(l) is 
a diagonal hermitian N x N periodic matrix, 

M(xi,x 2 ) =^V m (x 1 ,x 2 )r m (37) 

where functions rj m (xi, x 2 ) are real and periodic and will be determined through the remaining Bogomolnyi 
equations, equations (f3Tj) and (|32|) . Now, if we write Ai in (fl"8|) in terms of matrix M in the form 



At = iEij djMAr 1 , (38) 
Bogomolnyi equation (|34[) imposes to x the condition 

+ at t z ^ mTm X = , (39) 

where we have defined z — X\ + ix 2 . Remarkably, eq. (|39p is automatically satisfied by x as defined in 
eq. (|2"3")) . Hence, the only remaining task in order to have a complete solution to the Bogomolnyi equations 
is to determine M and P/. 



G 



5 Elementary U(l) x SU(2) non-Abelian vortex solutions in T 2 



In this section we complete the construction of string like solutions to the Bogomolnyi equations for the 
simplest U(l) x SU(2) gauge group case, leaving for the next section the extension to the general U(l) x 
SU(N) case. We shall consider the same number of flavors and colors, N = Nf — 2. We shall also restrict 
the analysis to the case of an elementary non-Abelian string solution with (n\,n 2 ) = (1,0). The case in 
which [n\,n 2 ) = (0, 1) is completely analogous. 

We have seen from the boundary conditions that the Higgs field $ can be factorized as a product of 
two functions, one (x) satisfying non-trivial boundary conditions, the other (A*), a strictly periodic function 
which remains to be computed. In the N = 2 case x is a 2 x 2 matrix satisfying (]22[) 

/ e i7r x| q\ 
X(xi+L 1 ,x 2 ) = I 1 jx(xi,x 2 ) 

fe^^T ()\ 

x(x\,x 2 + L 2 ) = I \x{x\,x 2 ) . (40) 



Then, we can write \ as 



X{xi) = (41) 

\ X22{X 1 ,X 2 )J 

with xn satisfying twisted boundary conditions 

Xn[xi+Li,x 2 ) = e L7 ' 1 ^Xii{xi-,x 2 ) 

Xu{xi,x 2 + Li) = e~™~xn{xi,x 2 ) (42) 

and X22 periodic in T 2 . 

These are very similar to the boundary conditions that arise in the Abelian Higgs model in the torus and 
this is the reason why construction of the solutions will closely follow [B]. Indeed, as in the Abelian case, we 
can find a solution to eqs. (|4*2"l) in the form 

i \ (■ (xi+ix 2 )x 2 \ 
Xn(xi,x 2 ) = exp I itt — I <d(x 1} x 2 ) , (43) 

with 

e(x 1 ,x 2 ) = 9 3 (, X -l±^ . (44) 

Here O3 is the Riemann theta function already defined in ()25|) . Then, x can he simply written as 



(fr (a ffi 3 )a3 )e(si,s2) o x 



X (x u x 2 )=^ XP V" ^ J " . (45) 

Without loosing generality we choose X22 = 1 and will accommodate A^ so that (f>f fulfills the appropriate 
boundary conditions. 

Concerning the gauge field, we shall take ansatz choosing M (see cq |57|) in the form 

M(xi,x 2 ) = exp {iriQTo + irjm) (46) 
With this choice one can see that the only non-zero components of the gauge field are A® and A\. 
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As explained at the end of the previous section, the ansatza? (p?Tj) and automatically solve one of the 
Bogomolnyi equations, namely eq. ([33|) . The problem is then reduced to solving the remaining two equations, 

F° = (,, f ,, f ■>),,, (47) 
F§ = i2«V/*Ve« . (48) 

where k is the ratio of the coupling constants, n — g/e. Here we have scaled coordinates and fields in the 
form 

Af^J^AT, tjS ^ y/ltf . (49) 
V V I 



Note that periods L\, L 2 are also rescaled according to 



L t ^J^-L t (50) 



Recalling equation (|36p . we write 

(j} i {x ll x 2 ) = X (xi,x 2 )M(x 1 ,x 2 )V f . (51) 
where the multiplets can be chosen in the form 

so that the two flavor components of the Higgs field becomes 

^,(*f'), ( 53 , 

with M11 and M22 the diagonal components of matrix M to be still determined. 

This, together with equation (flU)) . allows us to rewrite equations l[4"7|) - (|4"8j) in terms of strictly periodic 
fields 

ft = -(Af?\xn\ 2 M 2 n +^M 2 22 -2- 2 " 



F% = (^1x111^-^^ + ^-)^.. (54) 

Using (fl"9|) . ([38| and (|37|) . the periodic piece of the field strength can be written as 

F% = £i ,VV (55) 



where 



= £ij Vy , (56) 



V 2 = d 2 + 3 2 . (57) 

Then, Bogomolnyi equations (|5"4"|) become 

VV = -AA 2 |xii| 2 M 1 2 1 -AA 2 2 M 2 2 + 2-/f (58) 

VV = K 2 AA 2 |xii| 2 Mi 2 i-« 2 AA 2 M 2 2 + — - . (59) 
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or using equation (|57|) 

VV = -N^\ X ii\ 2 e~^ a ~^ -Nie-^°+^ + 2- (60) 

VV = K 2 AA 1 2 |xii| 2 e^(" - I ' 1 )- K 2 AA 2 2 e -(" +'" 1 ) + -^ . (61) 

L1L2 

From eqs. (|60l) -(|6T |) we see that 

(* - „) . wfia.iv>"-> + ^ (-L + 1) - 2 (62) 

If we now integrate both sides on torus, the l.h.s. vanishes since both 770 and rj\ are periodic. We then find 



J^dT (uf?\xu\ 2 e 
Now, since A/" 2 |xii| 2 exp(ryi — rjo) > one has 













- 


U 2 " 





(63) 



1 + 1 ) -1) / ,17, (04) 



LiL 2 \k 2 / / 7r 
or, calling A = L1L2 the area of the torus, 

A>*(l + ±>j=A c (65) 

That is, in order to have consistent solutions from our ansatz, there is a minimal critical area which we call 
j4 c , such that no solutions exists for A < A c . It will be convenient to introduce the parameter e, 

e = l-^ (66) 

which measures the departure from this critical area, with range < e < 1. 

We shall now solve the system (f6"0|) - (l6TT) and, for simplicity, we shall consider the case in which gauge 
coupling constants e and g coincides so that n 2 = 1 and the critical area reduces to A c = 2tt (the general 
case can be solved analogously). Defining 

h ± = ^( Vo ±m) (67) 

system ([oTJ|) - (|5T|) becomes 

V 2 h+ = 1-N%e~ 2h+ (68) 
V 2 /*- = e-Nl\ X ii\ 2 e- 2h " . (69) 

We shall construct vortex configurations starting from the trivial solution of eq. (|68l) 

Ml = 1 , h+ = . (70) 

Concerning eq.(|t)9"|). when e = 0, it also has a trivial solution 

Nf = , hT = 0. (71) 

which, together with (j70|) leads to 770=771 = 0. For £ ^ eq. ([69]) can be solved extending the method 
proposed in [6] for abelian vortices (see also [21], [22]) which consists in expanding h~ and Nf in powers of 
the Bradlow parameter e, and further expand the fields in fourier modes. The coefficients of the expansions 
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Figure 1: We plot F® 2 and 4>f4> for different Areas. When A = A c (e = 0) the solutions are trivial. When ^4^00 
(e — » 1) non-Abelian vortices on the plane are recovered. We plot solutions for different values of e. The area is 
written in units of -2^. 

obey recursive relations that allow the numerical calculation of the magnetic and Higgs fields. This is 
explained in an Appendix. 

Using this method, we have solved numerically eq. (|69p and obtained M\ and h~ for different values of e 
in the range < e < 1. From these results, the magnetic field F® 2 and the Higgs fields $ can be computed, 
using the equations 

F? 2 (xi,x 2 ) = V 2 h-( Xl ,x 2 ) (72) 

^rt-PM* 1 "^) j, (73) 

where 0(xi,X2) is defined in (|4~4")) . 

We show some of these solutions in Figure [T] In all these cases, no = 1 so the flux f T2 dTF^ = 27r . 
When varying the area A c < A < 00, solutions interpolate continuously from the trivial constant solution 
for e = to the non-Abelian vortices on the plane for e = 1. The vortex profiles are similar to those in the 
plane, with the magnetic field concentrated around the position of the vortex. At the center of the vortex, 
the upper component {<p l ) of the Higgs field, the one with a non-trivial winding, is zero, as it happens in the 
Abelian case. Typically, when the area is small A ~ A c the solutions converge fast, obtaining high precision 
by computing a few orders of the e expansion. In the infinite area limit A ^> A c the method converges much 
slower. In this case we have considered up to 40 orders of the expansion with more than 400 Fourier modes 
(which allows for a precision of less uncertainty than 10 -6 for the energy or magnetic flux). 

In the left panel of Figure [2] we show the magnetic field for the e — ► 1 case, which corresponds to an 
elementary ((1, 0)) non-Abelian vortex on the plane. The ansatz and the numerical method work as well for 
the study of multi vortex configurations, even when the vortices are not superimposed. We show in the right 
panel of Figure [5] a (2,0)-vortex configuration in the limit e — > 1. 
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Figure 2: We plot elementary and non-elementary 2-vortex configurations in the large area limit (^4 S> A c ). The 
model has gauge group (7(1) x SU(2). 

6 U(l) x SU(N) strings 

In this section we extend the analysis to the N > 2 case, with Nf — N. We start from the Bogomolnyi 
equations (|31[) - (j52")l and consider a non-elementary (n, 0, . . . , 0) vortex 

4 = - (tty - N) a, 

F™ = iVWK^t^dj , (74) 

where again k is the ratio of the coupling constants, k — g/e. Other elementary vortices like (0,n, . . . ,0), 
etc. can be analogously treated. 

We write the Nf = N flavor multiplets <j){ in the form 



<t>\ = Mexp('--l=(» ? -A/iV^l ?7 1 ))xn 



1 = --- = ^n = ^™p(-^=[ v » + ^= v l ) ) , (75) 



/27V V \fW^l 

where xn is defined in (|2"3")) . and 77 and r\ x (all other 77's are taken to be zero) are defined in (|57|) . With this 
ansatz, and defining 



we are able to rewrite equations (|T4"|) as 

V 2 »+ = (l-^^)-^exp(-^((iV-( 1 -^ + (l-K>-)) (77) 

\ L\L 2 N + >?L Y L 2 N J 
- Af?\xn\ 2 exp(-M((N-l)(l- K 2 )h+ + (l + K \N-l))h-y\ , (78) 



where 7 is defined in JTI 

As in the N = 2 case, we notice in equation ([75]) . that since 



dTV z h- = , (79) 

T 2 
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integrating equation ([75]) on the torus gives 



7° (jV-l)7° N 
LiL 2 N K 2 L t L 2 N ) J T 2 



1 + ^7+ by /I, ) / «n*>o, (80) 



and this again implies the existence of a critical area A c such that 

/r / N - 1 \ 

a ^ a ^^n{ 1 + — ) ' (81) 

We again consider the parameter e = 1 — A c /A in terms of which equations ([77|) - (|75|) read 

V2/l+ = (* + V- (!-" ' A/ " |CXP (-V^((^ ~ (! ~ " 2 )) fe+ + (! ~ ( 82 ) 

V 2 ft~ = e-AA 1 2 |xii| 2 exp(-VW((^-l)(l-^ 2 )^ + + (l + « 2 (^-l))^ _ )) • (83) 
For e = one has the simple solutions 

ft+ = ft - = . (85) 

Non-trivial solutions when e ^ can be obtained as before, Fourier expanding fields, and further expanding 

fields in powers of e. Order to order in e, one is left with recursive relations for the coefficients. These 
relations can be handled numerically as in the N = 2 case. 



7 Summary and discussion 

The main goal of this work was the study of field configurations corresponding to a periodic array of non- 
Abelian vortices. We have considered a Yang-Mills theory coupled to fundamental scalar matter, a model 
which can be seen as the truncated bosonic sector of a Af — 2 supersymmetric QCD. We have studied 
these configurations by solving the Bogomolnyi-Prasad-Sommerfeld equations of the theory. By analyzing 
the (twisted) non-trivial boundary conditions that the fields must satisfy on the two-torus, we were able to 
propose an ansatz that reduces the BPS equations to a a simpler set of ordinary non linear equations that 
can be solved numerically. These equations are solved perturbatively in powers of a parameter measuring 
the departure of the area of the torus from a critical minimal value. 

We have presented explicit solutions for the simplest gauge group U(2) which are the natural generaliza- 
tion of the ones studied by Gonzalez Arroyo and Ramos [5] to a non-Abelian Gauge theory. On the other 
hand, for large areas, our solutions converge to those studied in [TT]-[T3], for SUSY QCD. 

Our work could be extended in several directions. We have analyzed the case in which Nf — N . A 
natural extension would be to consider the case in which Nf > N to study non-Abelian semi-local strings 
[23] and this could be of interest in connection with low-energy effective actions for string theories. It is 
also natural to expect that the same ansatz presented here would work practically in the same way for 
Chern-Simons-Matter theories |24j-[25j. giving in this case origin to configurations of periodic, electrically 
charged non-Abelian vortices. Also, a similar analysis as the one presented here should be of use to study 
non-Abelian periodic vortex array configurations presenting BPS equations in the a non-Abelian model 
with adjoint matter [26], [28] or in the Standard Model [27]. The case considered here corresponds to the 
particular set of parameters dictated by supersymmetry and BPS equations. It is related, in the Abelian 
Higgs model, to the limit between Type I and Type II superconductivity, where vortices are non-interacting. 
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Away from this point, the full second order Euler Lagrange equations should be solved. This case, that 
would correspond to interacting vortices, is technically more involved to study. We expect that there exists 
a region in parameter space where the vortex- vortex interaction is repulsive giving rise to a lattice of vortices 
with a definite geometry. We hope to deal with some of this issues in the future. 
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Appendix 

To solve equation (JHU) 

V 2 h- =t-Nl\xn\ 2 e- 2h ~ , (86) 

we first define 

H(x 1 ,x 2 ) = e 2h ~^' x *\ (87) 
and rewrite (|86p in complex coordinates z = x\ + 1x2 

d z {H-%H) = l - (e-AfllxnfH- 1 ) . (88) 

Since for e = there is a trivial solution 

H = Constant, M = , (89) 

and considering that < e < 1, we can use s as a perturbative parameter and expand H and the normalization 
constant A/"i in powers of e 

00 00 00 

H = J2 H k s\ H ~ 1= Y. SkE ^ ^ = Y. Nk£k ■ ( 9 °) 

fc=0 k=0 k=0 

The coefficients Hk and -fffc are periodic functions and can then be Fourier expanded 

fj^ _ ^ ' ^(fe) e 27ri(nia:i/Li+n2X2/L2) _ ^ ' Jj(fe) e 27ri(rii3:i /£i+"22;2/£2) ^gjj 

niri 2 niri 2 

and the same can be done for |xn| 2 



IX11I 



= r\ n n e 27T i(nix 1 /L 1 +n 2 X2/L 2 ) ^g2) 



with normalized coefficients such that 7700 = 1. 

Inserting these expansions in eq. ([86|) one can determine order by order the coefficients, 



fr (0) = h {0) 

'"n 1 n2 "nin 2 



1 n\ — ri2 = 

n!^0,n 2 ^0 



where 



(1 ) I m = n 2 = 

"Tiin 2 " - nin 2 ' V 30 /* 



£ ni „ 2 = ttVt ( mi + — = J , (94) 



with T = L2/L1 the aspect ratio of the torus. In the same way one can calculate coefficients to any order Q 



m e 



ni = n 2 = 

KlL = { g w-c^-£g3 , (95) 

|g„ in j2 — «i / 0, n 2 ^ 
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with 

O-i 

9192 — ™l"2 " - 9l-ni,92-"2 ^9192 S?l-ni,?2-n a 

nin 2 fe=l 

Q-l 

Cflf*. = 277 E E ^Q-fc %x-n 1 , 92 -n 2 

niTi 2 fc=0 
Q-2 

9192 Z^ »i"2 'V— "1,92— "2 ^9192 Sgi-ni.ga-ns • 

Coefficients h nin2 , appearing in the expansion of H —1 , are obtained from the condition HH^ 1 = 1 

Q 

SCO) = V V /i (fc) h (Q ~ fc) (96) 

"9192 2-^1 2-^1 ™1™2 "9l-«1.92— "2 ■ \ JK) > 

nin.2 k=l 

One also has to find a recurrence relation for the coefficients Nf.- For this, the condition J T2 dTF 12 = 
implies 

N = 0, iV 1 = l 
O-i 

N Q = "EE V-n u -n 2 h[ k L N Q -k , Q > 1 ■ 
nin 2 fc=l 

(97) 

Computing these recursive relations we can obtain h~ and A/i, and with this compute the magnetic and 
Higgs field from eqs. ([72|) - (|751) . 
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